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We show that very light (50 - 90 GeV) axigluons with flavor- universal couplings of order g s /3 
may explain the anomalous top forward-backward asymmetry reported by both CDF and DO 
collaborations. The model is naturally consistent with the observed ti invariant mass distribution 
and evades bounds from light Higgs searches, LEP event shapes, and hadronic observables at the 
Z pole. Very light axigluons can appear as resonances in multijet events, but searches require 
sensitivity to masses below current limits. 
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1 Introduction 

The CDF and DO collaborations have recently re- 
ported measurements of the forward-backward asym- 
metry (Afb) in it production with intriguing devi- 
ations from the standard model prediction. CDF's 
result p] in the lepton plus jets channel reports an 
inclusive parton level asymmetry 



A FB (CDF) £j = (15.8 ± 7.4)% 



(1.1) 



If their measurement in the dilepton channel [2] is 
combined with this result, the asymmetry becomes 



A FB (CD¥)u+tj = (20.9 ± 6.6)% 



(1.2) 



and exceeds the standard model prediction ~ 5% 
[3]-[5] by more than 2 standard deviations. 

DO performs a similar search [6] in the lepton 
plus jets channel and reports an inclusive parton- 
level asymmetry 



A FB (D0)# = (19.6 ± 6.5)% 



(1.3) 



which is also more than 2a above the SM result. 
Taken together, these consistent deviations may be 
evidence for new physics in top quark production. 

While all the inclusive measurements are consis- 
tent with each other, CDF's lepton plus jets search 
sees sharp mass dependence [1] in the binned result 



A FB (M a < 450 GeV) 
A FB (M t t > 450 GeV) 



(-11.6 ± 14.6)% 
(47.5 ±11.4)% , 



where the high mass bin is 3.4 a above the SM pre- 
diction. Neither DO nor the complementary CDF 
dilepton search see the same effect; both find con- 
sistently positive > 2a deviations from the SM over 
the full M t i range. 

It has been observed that massive gluons with 
axial couplings can induce a large forward-backward 
asymmetry in ti production by interfering with stan- 
dard model processes [7]-jl9j. Motivated primarily 
by the mass dependent CDF result, these models 
predict asymmetries that rise uniformly with invari- 
ant mass and feature a sign flip near M t j ~ 450 GeV. 
Large (TeV scale) masses are typically required to 
satisfy dijet-resonance search bounds and suppress 
contributions to the ti invariant mass distribution. 
To produce an asymmetry with the observed sign, 
most models also require flavor violation and are 
severely constrained |20j by limits on flavor chang- 
ing neutral currents. For a comparison of heavy ax- 
igluons and other models that address the top asym- 
metry, see [21]. 

Relatively lighter axigluons (400 - 450 GeV) [22] 
can produce a large top asymmetry without flavor 
violation, but this mass scale is in tension with dijet 
resonance bounds and the differential M t i distribu- 
tion. Extra field content is generally required to 
broaden decay widths and avoid resonant enhance- 
ments to top quark observables. 

In this paper we propose a very light (50 - 90 
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Figure 1: Axigluon contribution to ti pair produc- 
tion. Interference with the standard model gluon 
exchange diagram generates Af nt . 

GeV), weakly coupled axigluon to explain the top 
asymmetry. The model inherits many of the fea- 
tures heavier axigluons enjoy, but counterintuitively 
avoids their experimental constraints by being light: 
dijet resonance searches suffer from large QCD back- 
grounds at low invariant masses, particles below the 
2mt threshold do not produce bumps in the ti invari- 
ant mass distribution, and nonresonant production 
suppresses new physics contributions to the ti cross 
section, which start at fourth order in the axigluon 
coupling. We find that the strongest upper bounds 
in this mass range come from Tevatron searches for 
light Higgs bosons produced in association with an 
additional 6-jet. The strongest lower bounds come 
from UA2 dijet searches and LEP measurements of 
the hadronic Z width. 

In Section [2] we describe our model; in Section 
[3] we discuss the details of our numerical simula- 
tion; in Section [4] we address the experimental con- 
straints; in Section [5] we compute the tt forward- 
backward asymmetry and compare theoretical pre- 
dictions with production- level data; in Section [6] we 
make some concluding remarks. 

2 Model Description 

We give the axigluon (G') flavor universal couplings 
to SM quarks 



CDg'G'^QT^Q 



(2.1) 



where g' = Xg s is the axigluon coupling constant, 
which we express in units of the strong coupling. 
This operator can arise from an extended SU(3)\ x 
SU(3)2 color group that breaks down to the diago- 
nal SU(3) C of QCD and gives rise to massive spin-1 
color octets |23j-|26j. For an axigluon of mass mc 
our effective model requires a UV completion at the 
scale 4vrm G //y = 1.7 TeV and 850 GeV for A = 0.3 
and 0.6 respectively. In this paper, we will focus 



only on the low energy effective theory and leave 
UV model building for future work. 

Without additional field content, all decays pro- 
ceed through operator in Eq. f|2. 1[) . so axigluons can 
only decay to quark pairs and give rise to dijet and 
four jet events for single and pair production, re- 
spectively. Since we work in the regime where the 
axigluon is below the ti threshold, the total width 
is [271 



rif 2 
r G / = -r-a s X m G > 
b 



(2.2) 



where nj is the number of active fermion flavors. 
For m G i = 80 GeV and A = 0.4, this width is Y G > ~ 
1.1 GeV. 

The differential cross section for the process qq — > 
ti in the CM frame is a sum of standard model, in- 
terference, and axigluon terms 



where 
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Here f3 = \/l — 4to 2 / s is the top quark velocity and 
9 is the angle between the incoming quark and out- 
going top in the CM frame. A forward-backward 
asymmetry can only arise from terms with odd pow- 
ers of cos#, so the effect is due entirely to inter- 
ference. In the presence of both vector and axial- 
vector couplings, there is an additional small con- 
tribution to the asymmetry from the new-physics 
squared term. 

Note that the asymmetry generating term Af nt 
is proportional to (s — m G ,). For heavier axigluons, 
this dependence gives rise to a negative asymme- 
try because the mass is typically larger than the 
partonic CM energy. To compensate, many models 
introduce opposite sign couplings to the first and 
third generations. In our case, tuqi < s for on-shell 
ti production, so the asymmetry is always positive 
and flavor violation is unnecessary. 
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3 Simulation and Acceptances 

In the lepton plus jets analysis, CDF unfolds raw 
data by deconvolving their detector simulation and 
jet algorithm to yield a partonic data set from events 
that survive cuts at the detector level. To compare 
our model predictions with this data, it is neces- 
sary to generate an event sample with partonic ti 
pairs in the final state. However, knowing the pre- 
dicted cross section and experimental luminosity is 
not enough to properly normalize kinematic distri- 
butions from the partonic simulation; we must also 
know the detector level acceptances. We thus per- 
form two simulations: one at the partonic level to 
make our plots and one at the detector level with 
CDF's cuts to compute the acceptances that nor- 
malize these distributions. 

We simulate the partonic process pp — > ti in 
MadGraph 5 [29j using a model file generated with 
FeynRules [3D]. This file adds the operator in Eq. 
(|2.ip to the full standard model Lagrangian so that 
the process in Figure [1] contributes to ti produc- 
tion and gives rise to interference with SM gluon- 
exchange. 

For the acceptances, we also perform a more re- 
alistic simulation (pp — > ti — > iv + Aj) using Pythia 
|31| for the parton shower and PGS [32J for detec- 
tor effects. To compare with CDF's lepton plus jets 
search, we impose the following cuts: at least four 
jets with Et > 20GeV and at least one 6-tag; for 
non-6 jets \rjj\ < 2, for 6-jets |%j| < 1; large missing 
energy Et > 20GeV; and exactly one electron or 
muon with p T > 20GeV and |%| < 1. 

Note that there is some error introduced by this 
approximate method. A complete comparison with 
experimental data would not only run a full detec- 
tor simulation (including Pythia and PGS), but 
also identify top quarks with a least-squares kine- 
matic fit and unfold the detector-level output us- 
ing the CDF algorithm that reconstructs partonic 
events from raw data. Nonetheless, our approach 
accurately reproduces CDF's standard model expec- 
tation for the ti invariant mass distribution^ so the 
error introduced by a constant acceptance function 

1 Although the forward-backward asymmetry arises only 
at loop level in the SM, its numerical value is tiny (~ 5%), so 
this tree level method also adequately reproduces the (nearly 
symmetric) SM predictions for the Ay = yt — Vi rapidity 
distributions in pQ. 



is likely to be small in our case as well. We leave 
the full unfolding for future work. 

4 Experimental Constraints 

Models that explain the top asymmetry must agree 
with the ti invariant mass distribution and total 
cross section, both of which are in good agreement 
with standard model predictions. Any candidate 
model with an s-channel mediator must satisfy con- 
straints from dijet resonance searches at hadron col- 
liders. In our case, we must also contend with a va- 
riety of older measurements that set lower bounds 
on new colored particles. 

4.1 Top Quark Measurements 

The ti cross section at the Tevatron has been mea- 
sured to be o"jf P ' = 7.50 ±0.48 pb [33J , which agrees 
with the standard model prediction in perturbative 
QCfH, ~ (6.32 - 7.99) pb for m t = 172 GeV 
[34"] . The leading order result, (a^f 1 )lo — 5.63 pb, 
computed with MadGraph, implies a SM i^-factor 
between 1.12 and 1.42. 

Including an axigluon with mg' = 80 GeV and 
A = 0.4, gives a total LO cross section of (cS^lo = 
6.08 pb, which is only an 8% increase over the SM 
LO result. This minor enhancement is due entirely 
to Aa xi in Eq. (|2.ip , which is fourth order in the ax- 
igluon coupling; the interference term Af nt does not 
contribute to the total cross section. Although com- 
puting higher order corrections is beyond the scope 
of this work, the color structure of the axigluon ex- 
change diagrams is identical to that of the relevant 
SM processes, so we expect higher order corrections 
to be of similar magnitude, though a more precise 
calculation is necessary to take into account the ad- 
ditional interference. As long as the K factor does 
not differ substantially from that of SM production, 
the total ti cross section stays in good agreement 
with experiment. For the remainder of this paper, 
we will assume the K factor to be 1.2, so our bench- 
mark cross section becomes 7.3 pb. 

For very light axigluons (mc <C 2mt), top pair 
production is nonresonant, so the invariant mass dis- 
tribution is also in good agreement with experiment. 
In Figure [2] we show the simulated M t i distribution 

2 For complementary calculations see [351136] . 
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Figure 2: Tevatron invariant mass distribution for 
tt pairs (blue, color online) including both axigluon 
and background contributions. Data points and 
standard model background (purple) are taken from 
CDF's lepton plus jets search pQ. Here we use 
A = 0.4 and mc = 80 GeV. After including a a 
if-factor of 1.2, the top cross section is a t i = 7.3 
pb. Applying the CDF cuts (see Section [3]) gives an 
acceptance of 2.6%. 



(blue) plotted alongside the CDF data points and 
standard model background (purple) taken from the 
lepton plus jets search pQ. 

4.2 Dijet Resonance Searches 

Quark coupled axigluons give rise to two and four jet 
events from single and pair production, respectively. 
Our mass range of interest (50 - 90 GeV) is safe from 
Tevatron (33 EH] and LHC gl, [40] dijet resonance 
searches, which do not set bounds on masses below 
180 and 200 GeV, respectively. A preliminary AT- 
LAS analysis of multijet events [UJ sets limits on 
color octet scalars with narrow widths, but does not 
constraint masses below 100 GeV. With lower search 
thresholds, this model may be testable at both the 
Tevatron and LHC, however, signal and background 
are expected to be large at both colliders [4*2] . 

The UA2 search for hadronic W and Z decays 
[4*3] measures the exclusive two-jet mass spectrum 
between 48 and 300 GeV, which constrains the light 
axigluon parameter space. Using 4.7 pb" 1 for Mjj > 
66 GeV (and 0.58 pb" 1 for 48 GeV < Mjj < 66 
GeV), the combined W and Z resonances are ex- 
tracted with a bi-gaussian fit above a smooth back- 
ground function normalized to the data. The best fit 



bi-gaussian signal spans the Mjj range between 70 
and 100 GeV and yields a cross section of a ■ Br(W, Z 
jj)obs. = 9.6 ± 2.3 ± 1.1 nb, whose central value 
exceeds the SM prediction at NLO, a ■ Br(W, Z —> 
jj)sM = 5.8 nb, by almost a factor of two. 

Although a three-gaussian fit and a QCD back- 
ground prediction are necessary to properly con- 
strain axigluons using this data, we can extract a 
rough bound by finding (\,mc>) values for which 
the combined SM and new-physics predictions ex- 
ceed the observed number of events under the best 
fit gaussian by 2a. In Figure [3] we plot the exclu- 
sion boundary (yellow dot-dashed line) determined 
using Madgraph, Pythia, and PGS to simulate our 
signal. 

For dijet masses below 70 GeV, the UA2 analy- 
sis does not attempt to fit any signal, so a possible 
resonance would almost certainly have been missed 
given the very large background in this mass range. 
Even near mw and m?, the signal/background ratio 
is only a few percent and the gauge boson peak is not 
visible to the naked eye (see Figure 5 in [33]) prior 
to a rescaling that emphasizes the region around 
the known W and Z masses. Since the background 
model for this search is purely data-driven, the low- 
mass region does not impose a meaningful constraint 
without a dedicated bump hunt. 

4.3 Light Higgs Searches 

Tevatron searches that look for light Higgs bosons 
produced in association with 6-jets {pp — > hb — > bbb) 
are sensitive to axigluon decays into 6-quarks. Since 
these searches require at least three 6-tags to reduce 
the QCD multijet background, the bounds they im- 
pose on a(hb) • Br(h — > bb) also apply to the pro- 
cesses pp — ^ G'b — ^ bbb and pp — > G'bb — >■ bbbb, the 
latter of which can also arise from pair produced ax- 
igluons. However, the CDF [H] and DO [25] results 
only apply to masses above 90 GeV; light axigluons 
fall below the sensitivity threshold. To be conser- 
vative, we will only consider masses below 90 GeV 
where the 36 constraints do not apply. 

The authors in |46j use Tevatron Higgs searches 
in the associated production channel, pp — > Wh — > 
(lv){bb) to exclude axigluons with A = 1 between 
75 — 125 GeV assuming Br(G' — > bb) = 1/5. In our 
case with A = 0.4, the Tevatron qq — > WG' cross 
section decreases by a factor of A 2 , which reduces 
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Figure 3: Allowed axigluon parameter space in 
the (A, mc) plane plotted alongside bounds from 
dijet-resonance searches and T(Z — > hadrons) mea- 
surements assuming different extractions of a s . The 
blue and purple bands (color online) are regions 
favored by the combined CDF/DO inclusive asym- 
metry measurements at la and 2a, respectively. 
The dot-dashed yellow curve marks the approxi- 
mate 2a bound above which model predictions ex- 
ceed UA2 dijet limits from hadronic W and Z decays 
(see Section 14. 2j) . The solid black curve marks the 
boundary above which corrections to the hadronic 
Z width exceed the observed value by 2a assuming 
the standard model extraction of a s (mz) = 0.1184. 
The dashed and dotted black curves give the same 
bound, but respectively assume 2.5% and 5% re- 
ductions to the SM value of a s (mz)- Reductions 
of this magnitude are typical of light axigluon con- 
tributions to the QCD beta function (for a discus- 
sion see Sections 14.51 and 14. 6j) . The region above 
rriQi > 90 GeV is excluded by Tevatron 36-searches. 
Since LEP event shapes rule out gluon-coupled ad- 
joint fermions around 50 GeV, our model may en- 
counter a stronger lower bound since axigluons also 
couple to quarks, but a proper analysis is necessary 
to set the correct limit. 

the axigluon signal a ■ Br from ~ 50 pb down to ~ 5 
pb for rriQ> = 50 GeV also assuming Br(G' — > bb) = 
1/5. This falls safely below the quoted bound of 
< 20 pb, however, this number is based on analysis 
from an unpublished talk, so its status is not clear. 
Current Tevatron searches for the associated pro- 
duction of Higgs bosons are not sensitive to masses 
below 100 GeV [471 [48] . 

Naively it would appear that LEP searches in the 



Higgstrahlung channel [49]- [52] e + e~ — > Zh — > 4j 
would be sensitive to light axigluons produced in 
e + e~ — > Z* — > qqG' —> 4j events. However, the 
event selection algorithms in these analyses look for 
kinematics that fit the Higgstrahlung topology in 
which the invariant masses of jet pairs produce both 
Z and Higgs resonances. In events with on-shell 
axigluons, all four jets arise from virtual Z exchange, 
so this possibility is highly disfavored. Furthermore, 
this process occurs at order A 2 and suffers additional 
phase-space suppression. 

Similar considerations apply to LEP measure- 
ments of triple gauge boson couplings [53j-[56j which 
look for e + e~ —> W + W~ , ZZ — > 4j events. These 
analyses select events using neural network algo- 
rithms designed to identify diboson production; light 
axigluons arising from Z exchange have very dif- 
ferent kinematics and fail this selection, which re- 
quires some combination of jet pairs to reconstruct 
at least one gauge boson mass. At the higher end 
of our mass range (mc > 80 GeV) it may be possi- 
ble for an axigluon to fake a hadronically decaying 
SM gauge boson, but the other two jets would not 
reconstruct a resonance. The coupling and phase- 
space suppression also diminish the rate at these 
searches, so axigluon production is negligible com- 
pared to tree-level diboson and QCD background 
processes. 

4.4 Event Shapes 

Constraints on light colored-particles have been ex- 
tracted from the analysis of event shapes at LEP. 
Comparing multijet data with calculations in soft 
colinear effective theory (SCET) rules out color ad- 
joint fermions below 51 GeV at 95% confidence |57j . 
However, this approach assumes that the new field 
couples only to gluons, with no tree-level quark in- 
teractions. To set a proper lower bound, it is nec- 
essary to repeat this analysis with more general as- 
sumptions, however, it is unlikely that this would 
yield a more lenient limit so we will not consider 
masses below 50 GeV. 

LEP studies of four-jet events from Z decays 
[SH]-[SI] can be sensitive to light, colored particles 
that couple to quarks. Various angular distributions 
are used to successfully distinguish SU(3) C QCD 
from alternative abelian theories of the strong force, 
so the presence of light axigluons could potentially 
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spoil this success. However, using Madgraph to gen- 
erate four-jet Z decays at the parton level, we find 
that the presence of an axigluon (A = 0.4) in our 
mass range does not qualitatively distort these an- 
gular distributions relative to the QCD prediction. 
This is unsurprising since 0(10%) of SM hadronic 
Z decays produce four-jets - the exact number de- 
pends on y cut and other jet algorithm details [62] 
- whereas in our model only 0(0.1%) of hadronic 
decays proceed through Z — > qq G' — > 4j prior to 
imposing cuts (see Section I4.6p . For higher ener- 
gies probed by LEP II (y/s « 200 GeV), the total 
e + e~ — > Z* — > qqG' — > 4j rate is similarly negligible 
compared to SM four jet production; this conclusion 
is robust for values of y cu t spanning several orders 
of magnitude. 

4.5 Running of a s 

Since axigluons couple to the strong sector, they give 
rise to loop diagrams that modify the QCD beta 
function above the scale mc ■ The standard model 
running between energy scales Q and [x is given by 



a s (Q 2 



ot s (l^ 



l + fea s ( M 2 )log(g 



(4.1) 



where, to leading order, b = (33 — 2rif)/12-K and 
nj is the number of active flavors. Since axigluons 
have the same quantum numbers and self couplings 
as gluons, their principal effect on the running is to 
double the gluon contribution to the beta function 
above vtiqi: 6 — >• (2x33— 2n/)/127r. This accelerates 
asymptotic freedom and yields smaller values of a s 
near the weak scale. 

While this adjustment naively jeopardizes the 
agreement between theory and experiment for the 
running, the experimental extraction of a s depends 
entirely on the assumed validity of standard model 
QCD with no additional field content [63J. At each 
energy scale, an a s -dependent observable is equated 
to the SM prediction and the resulting data point is 
extracted implicitly. If light new states were present 
in the strong sector, this data would completely ig- 
nore their contributions, so the current agreement 
between theory and experiment does not constrain 
our model. 

To roughly estimate the axigluon correction to 
Oi s (jnz), we use a well-measured value of a s below 



rriQi as an IR boundary condition and evolve it with 
the new beta function. This method is crude be- 
cause even low-energy observables used to extract 
a s depend somewhat on virtual axigluon processes, 
which are ignored in the extraction of reported mea- 
surements. Nonetheless, using the boundary con- 
dition a s (14.9GeV) = 0.160, [63] the weak-scale 
value becomes a s (mz) = 0.105,0.110, and 0.115 for 
mc = 50, 65 and 80 GeV, respectively. Different 
IR boundary conditions give similar downward cor- 
rections of order a few percent relative to the SM 
extraction a s (mz) = 0.1184. Note that this result 
is independent of A since axigluons couple to gluons 
with QCD strength. 

This model also predicts a kink in the running 
of a s near ma> ■ Our mass range of interest (50 - 90 
GeV), however, overlaps with a region where data 
points are sparsely distributed with relatively large 
error bars (see Figure 6 in |63J) compared to the data 
set as a whole. Kinks in the slope of a s would, there- 
fore, be unlikely to stand out in the data. Nonethe- 
less, a model-dependent extraction of a s is neces- 
sary to evaluate the possibility of kinks or overall 
data shifts due to new physics contributions. 

4.6 Hadronic Z Width 

The strongest lower bound on mp comes from vir- 
tual and three-body corrections to the hadronic Z 
width. Axigluons that couple to quarks with QCD 
strength (A = 1) enhance this width by a factor of 



1 + — / {m z /m G ,) + 0(a\ 

IT 



(4.2) 



where / is a function derived in [631 [65] . The LEP 
measurement of T(Z — > hadrons) and the extracted 
value of a s (mz) constrain the size of f (mz/mG') 
and severely restrict axigluon masses: ma 1 > 570 
(365) GeV for A = 1 at the 65% (95%) confidence 
level 06]. 

However, / is highly nonlinear, so the mass con- 
straint is extremely sensitive to the axigluon cou- 
pling. In our scenario, the constraint on / applies 
to the combination A 2 /, which dramatically weakens 
the lower bound on mc- Furthermore, following the 
discussion in Section [4"3l light axigluon (ma 1 < mz) 
contributions to the QCD beta function generically 
decrease the value of a s (mz) at the percent level. 
Since this is used to compute QCD corrections to the 
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SM prediction for T{Z — > hadrons) [66| . a smaller 
value opens up more allowed parameter space for 
new physics; the positive axigluon contribution to 
the width compensates for a slightly smaller SM re- 
sult which is reduced by the new value of a s . 

In Figure [3] we plot 2a exclusion bounds from 
the hadronic Z width on the (A,mc) plane along- 
side the regions favored by combined CDF and DO 
Afb measurements (discussed in Section [5]) . The 
solid black curve uses the standard model extrac- 
tion a s (mz) = 0.1184±0.0007[63] and the measured 
T(Z hadrons) = 1.744+0.002 GeV [67] to identify 
parameters for which the theoretical prediction ex- 
ceeds the measured central value by 2a. Also plotted 
are the 2a bounds assuming 2.5 % (black dashed) 
and 5% (black dotted) reductions in a s {mz) due to 
the modified running that includes axigluon contri- 
butions. These curves show how sensitive the bound 
is to modifications in A and a s (mz)- Since we gener- 
ically expect light axigluons to reduce the value of 
oisijnz) by a few percent relative to the SM extrac- 
tion, the dashed and dotted curves are more faithful 
to the underlying physics. Given the sensitivity of 
the bound, a proper extraction of a s involving ax- 
igluon processes is necessary to accurately constrain 
the parameter space; the limits in Figure [3] serve 
merely to illustrate the impact on the allowed re- 
gion. 

4.7 Bounds from a(e + e~ — > hadrons) 

The authors in [63] calculate H axigluon corrections 
to the ratio 



R(s) 



a(e + e — > hadrons) 
e 4 /12vrs 



(4.3) 



at the scale ^/s = 34 GeV and thereby exclude 
masses below 50 GeV at 95% confidence assuming 
A = 1. As with the hadronic Z width, the correc- 
tions for this process are proportional to the factor 



in Eq. (j4.2j) with the replacement a s — > A a s , so 
the discussion in Section 14.61 applies to this bound 
as well. Since T(Z — > hadrons) is extracted from 
R data at the Z pole, the allowed parameter space 
in Figure [3] is automatically consistent with bounds 
from R near yfs = mz- For smaller energies in our 



3 Note that [64| corrects some minor, yet consequential er- 
rors from an earlier paper [5S] that placed a far stronger lower- 
bound on the mass. 
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Figure 4: Inclusive top anti-top rapidity difference 
distribution plotted against unfolded CDF data. 
Here we use the same model parameters as in Fig. [21 
The blue histograms include both signal and stan- 
dard model background. Both data and background 
(purple) are taken from pQ. This plot omits the 
small, loop level asymmetry generated by SM pro- 
cesses. 

range of interest, yfs G 50 — 90 GeV, the uncertain- 
ties on the R data are larger than those at the Z 
pole [C>7| . so the bound is weaker. 



5 Forward Backward Asymmetry 



The forward-backward asymmetry can be written 



A 



FB 



N(Ay > 0) - N(Ay < 0) 
N(Ay > 0) + N(Ay < 0) 



(5.1) 



where Ay = yt—yi is the rapidity difference between 
the top and anti-top quarks. 

In Figure [3] we show the favored parameter space 
in the (A, rriQi) plane. The blue (purple) band repre- 
sents the region of la {2a) agreement with the com- 
bined CDF, Eq. CCD, and DO, Eq.CCI inclusive 
measurements. For typical points in these regions, 
the model predicts a positive asymmetry of order 
20%. 

In Figure U] we show the inclusive ti rapidity- 
difference distribution plotted against the CDF data. 
The signal simulation is identical to that used to 
generate Figure [2] with txiqi = 80 GeV and A = 0.4. 
After applying the cuts described in Section [31 the 
acceptance is is 2.6%. This plot only depicts the ef- 
fects of tree- level processes; the histograms do not 
include the small asymmetry induced by standard 
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model processes. However, the numerical results in 
Fig. [3] include the full asymmetry with both SM and 
new physics contributions. 

Although our simulation gives an acceptable fit 
to the rapidity data, some of the bins are more than 
la away from data points. We, however, do not 
expect perfect agreement at this level of analysis. 
The distribution in Figure H] is a rough approxima- 
tion of the full theory prediction which requires both 
a full CDF detector simulation and the subsequent 
unfolding for a proper comparison with data. 

In Figure [5] we show the theory prediction for 
the mass dependent asymmetry ApB{M t i) plotted 
alongside the unfolded CDF data. Like other light s 
channel mediators, light axigluons predict a positive 
asymmetry throughout the whole range of invariant 
masses. While the agreement at low invariant mass 
is not ideal, neither DO nor the CDF dilepton mea- 
surement observe strong mass dependence, so the 
significance of the mass-dependent data is not clear. 

Note that in Figures [2j |4] and [5] we only com- 
pare the model to CDF results because their pub- 
lished distributions feature production-level data, 
which allow for a direct comparison with parton 
level simulations. Comparison with DO's distribu- 
tions requires a detailed understanding of their de- 
tector simulation, which is beyond the scope of this 
work. Our conclusions have emphasized inclusive 
results from both collaborations since these are in 
better agreement with each other than the more con- 
troversial mass-dependent data. 

6 Conclusions 

We have shown that a light axigluon with flavor 
universal couplings can generate a large, positive tt 
asymmetry and naturally agrees with measurements 
of da/dM t t. The model has viable parameter space 
consistent with light Higgs bounds, dijet resonance 
searches and measurements of the hadronic Z width. 

For masses between 50 — 90 GeV and quark cou- 
plings in the range 0.3 g s — 0.6 g s , the theoretical pre- 
diction for the parton-level top asymmetry is in good 
agreement with inclusive results from both CDF and 
DO. The asymmetry is proportional to (s — niQ,), so 
the sign of Afb is always positive for on shell top 
pair production with y/s > 2m t 3> mg< . 

In the presence of a light axigluon, both the pre- 





i 1 


» 









400 500 600 700 BOO 



Figure 5: Theory prediction for the mass dependent 
ti asymmetry (purple histograms) plotted against 
the binned, unfolded CDF data in the lepton plus 
jets channel [T]. Here we use the same model pa- 
rameters as in Fig. [2 For comparison with CDF, 
the bin sizes are 50 GeV for M a < 600 GeV and 100 
GeV for larger invariant masses. Since the interfer- 
ence term in the differential cross section, Eq. (|2,5p . 
is proportional to (s — itIqi), the asymmetry is al- 
ways positive for on-shell tt production. This is a 
generic feature of light axigluon models. 

dieted and observed values of a s are modified at the 
percent level. A reanalysis of a s (y/s) measurements 
could reveal small downward shifts in the data since 
the modified beta function accelerates the running 
of a s in the presence of an axigluon. The down- 
ward shift in a s also decreases the SM predictions for 
T(Z —7- hadrons) and a{e + e~ —¥ hadrons), which ex- 
pands the parameter space for (A, me) values that 
explain the top asymmetry. 

Although the QCD background at low masses is 
formidable, it may be possible to revisit UA2 dijet 
data and perform a dedicated bump hunt in the low 
mass region with updated background calculations. 
It should also be possible to include light axigluons 
in a SCET reanalysis of event shapes in LEP data, 
which would likely set the strongest lower bound on 
this model. 

If very light axigluons explain the top forward- 
backward asymmetry, the Tevatron and LHC ex- 
periments should, in principle, be able to observe 
resonances in two and four jet events from single 
and pair production. Since the effective model pre- 
sented in this paper demands a UV completion at 
energy scales near the LHC's designed sensitivity, 
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we predict new physics around the TeV scale, but 
the specific signals are model dependent at this level 
of description and would be interesting to pursue as 
future work. 
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